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Recently, nuclear vibrational contribution signatures in 2D electronic spectroscopy 
have attracted considerable interest, in particular as regards interpretation of the 
oscillatory transients observed in light-harvesting complexes. These transients have 
dephasing times that persist for much longer than theoretically predicted electronic 
coherence lifetime. As a plausible explanation for this long-lived spectral beating in 
2D electronic spectra, quantum-mechanically mixed electronic and vibrational states 
(vibronic excitons) were proposed by Christensson et al. [J. Phys. Chem. B 116, 
7449 (2012)] and have since been explored. In this work, we address a dimer which 
produces little beating of electronic origin in the absence of vibronic contributions, 
and examine the impact of protein-induced fluctuations upon electronic-vibrational 
quantum mixtures by calculating the electronic energy transfer dynamics and 2D 
electronic spectra in a numerically accurate manner. It is found that, at cryogenic 
temperatures, the electronic-vibrational quantum mixtures are rather robust, even 
under the influence of the fluctuations and despite the small Huang-Rhys factors of 
the Franck-Condon active vibrational modes. This results in long-lasting beating 
behavior of vibrational origin in the 2D electronic spectra. At physiological temper¬ 
atures, however, the fluctuations eradicate the mixing and, hence, the beating in the 
2D spectra disappears. Further, it is demonstrated that such electronic-vibrational 
quantum mixtures do not necessarily play a signihcant role in electronic energy trans¬ 
fer dynamics, despite contributing to the enhancement of long-lived quantum beating 
in 2D electronic spectra, contrary to speculations in recent publications. 
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I. INTRODUCTION 


Recent observations of long-lived beating phenomena in the two-dimensional (2D) elec¬ 
tronic spectra?^ of photosynthetic pigment-protein complexes have stimulated a consider¬ 
able increase in activity in the interdisciplinary field of molecular science and quantum 
physics.^I®^ It has generally been assumed that electronic coherence decays sufficiently 
rapidly that it does not affect the nature of electronic energy transfer in photosynthetic 
light harvesting systems. Nevertheless, these recent experiments have attributed long-lived 
quantum beats to an electronic origin. Initially, 2D electronic spectroscopic experiments 
were conducted for Fenna-Mathews-Olson (FMO) pigment-protein complexes isolated from 
green sulfur bacteria at a cryogenic temperature, 77 and revealed the presence of 

quantum beats persisting for a minimum of 660 fs.^ However, it is generally thought that 
coherence at physiological temperatures is fragile compared to that at cryogenic tempera¬ 
tures, because environmental fluctuation amplitudes increase with increasing temperature. 
To clarify this issue, electronic coherence lifetimes in the FMO complex were theoretically 
examined, and it was predicted that electronic coherence in this complex could persist for 
700 and 300 fs at 77 and 300 K, respectively.^ These theoretical predictions were consistent 
with newer experimental results at physiological temperatures for the FMO complex.® How¬ 
ever, Panitchayangkoon et alW demonstrated that quantum beats in this complex persist 
for 1.5 ps at cryogenic temperatures, although the theoretical model does not produce elec¬ 
tronic coherence with this lifetime. Hence, nuclear vibrational contribution signatures in 
2D spectroscopy have attracted considerable interest of late, in particular as regards inter¬ 
pretations of the oscillatory transients observed in light-harvesting complexes, which persist 
for significantly longer than the predicted electronic dephasing times. Several independent 
publications have alluded to nuclear vibrational effects as an explanation for this long-lived 
beating behavior.^^^®^ 

In photosynthetic pigment-protein complexes, the Huang-Rhys factors of chlorophyll 
(Chi) and bacteriochlorophyll (BChl) are generally thought to be small,®®^ suggesting that 
the photophysics therein is primarily electronic in nature, rather than vibrational. Indeed, 
the recent 2D electronic experiments on BChl molecules in solution did not hnd signihcant 
vibrational coherence.^ In relation to this, Christensson et proposed that resonance be¬ 
tween electronic states and active Franck-Condon vibrational states serves to create vibronic 
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excitons, i.e., quantum mechanically mixed electronic and vibrational states. Such states 
have vibrational characteristics and enhanced transition dipole moments owing to inten¬ 
sity borrowing from the strong electronic transitions. Hence, it was suggested that coherent 
excitation of the vibronic excitons produces oscillations in the 2D signal that exhibit picosec¬ 
ond dephasing times. Tiwari, Peters, and Jona^ commented that the excitonically mixed 
electronic and vibrational states also lead to an enhancement of the vibrational coherence 
excitation in the electronic ground state. It was also argued that this effect could explain the 
long-lived oscillations in the FMO complex. Regarding electronic energy transfer dynamics, 
Womick and MorarpSI demonstrated that a vibrational-electronic resonance enhances the 
energy transfer rate in a cyanobacterial light-harvesting protein, allophycocyanin, with the 
use of the vibronic exciton model. 

Quantum mechanically mixed electronic and vibrational states, or vibronic excitons, are 
plausible as an explanation for long-lived spectral beating in 2D electronic spectra. How¬ 
ever, a question naturally arises concerning the interplay between the electronic-vibrational 
resonance and electronic energy fluctuations induced by the environmental dynamics.^^ In 
general, the energy eigenstates due to the quantum mixing of electronic and vibronic exci¬ 
tations are obtained via diagonalization of the Hamiltonian comprising the Franck-Condon 
transition energies and electronic interactions. It should be noted that these are indepen¬ 
dent of any environmental factors such as temperature, reorganization energy, and fluctua¬ 
tions. Ishizaki and Flemin^^^IlZl characterized the impact of environmental factors upon the 
quantum delocalization using the concurrence .^^^21 They visually demonstrated that smaller 
electronic coupling, larger reorganization energy, and higher temperature cause the dynamic 
localization,®^ even in the case in which two electronic states resonate in a coupled homo¬ 
dimer. Recently, IshizakP^ explored the influence of environmentally induced fluctuation 
timescales on the quantum mixing between electron donor and acceptor states in a model 
photo-induced electron transfer reaction. It was demonstrated that fast fluctuation and 
correspondingly fast solvation eradicates the quantum mixing between the donor and ac¬ 
ceptor in the vicinity of the crossing point of the diabatic free energy surfaces, leading to 
a Marcus-type nonadiabatic reaction, whereas slow fluctuation sustains the quantum mix¬ 
ing and prompts the electron transfer reaction in an adiabatic fashion.®' Therefore, it is 
natural to question whether the nature of the quantum mechanically mixed electronic and 
vibrational states is altered through dynamic interactions with the environment. 
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The main purpose of this paper is to explore the impact of environmentally induced 
fluctuations upon the quantum mechanically mixed electronic and vibrational states in a 
simple electronic energy transfer system through electronic energy transfer dynamics and 2D 
electronic spectra calculations. The electronic energy transfer system used here is a coupled 
hetero-dimer, modeled on BChl 3 and 4 in the FMO protein. Further, we investigate the 
extent to which vibrational modes play a role in electronic energy transfer dynamics under 
the influence of the environmentally induced fluctuations. 


II. MODEL 


To clarify the impact of the protein-induced dynamic fluctuations upon quantum mechan¬ 
ically mixed electronic-electronic and electronic-vibrational states in a systematic fashion, we 
consider the simplest electronic energy transfer system, a coupled hetero-dimer. To describe 
electronic energy transfer (EFT), we restrict the electronic spectra of the m-th pigment in 
a pigment-protein complex (PPG) to the ground state, and the hrst excited state, 

|(Pme), although higher excited states are sometimes of consequence in nonlinear spectro¬ 
scopic signals.^22^ Thus, the model Hamiltonian of a pigment-protein complex comprising 
two pigments is expressed as 

2 2 2 

-^PPC = EE ma ){^ma\ + EE ^Jmnl'^me) {'^mg\ ® ng) ne\■ 

m=l a=g,e m=l n=l 

Here, Hmai^m) represents the diabatic Hamiltonian for the environmental and nuclear de¬ 
grees of freedom (DOFs), Xm, when the system is in the \(pma) state for a = g,e. The 
electronic coupling between the pigments, hJmn, may also be modulated by the environmen¬ 
tal and nuclear DOFs. In the following, however, we assume that the nuclear dependence 
of hJmn is vanishingly small and employ the Condon-like approximation as usual. 

The Franck-Condon transition energy of the m-th pigment is obtained as 




mg-, 


( 2 . 2 ) 


where the canonical average has been introduced, (... )ma = Tr(... /5^„) with the canonical 
distribution, = e“^^"*“/Tr Here, /3 is the inverse temperature, 1/k-QT. The 

electronic energy of each diabatic state experiences fluctuations caused by the environmental 
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and nuclear dynamics; these dynamics are described by the collective energy gap coordinate, 
such that 


By dehnition, the coordinate, itm, includes information associated not only with the elec¬ 
tronic excited state, but also the electronic ground state. In this work, we assume that the 
environmentally induced fluctuations can be described as Gaussian processes and that the 
relevant nuclear dynamics are harmonic vibrations. Under this assumption, the dynamic 
properties of the environmental and intramolecular vibrations are characterized by several 
types of two-body correlation functions of Um{t) = as shown below. 

The environmental dynamics and intramolecular vibrational motion affecting the elec¬ 
tronic transitions can be characterized by the nonequilibrium energy difference between the 
electronic ground and excited states as a function of a delay time, t, after the photoexcitation.^ 
The linear response theory allows one to express the nonequilibrium energy difference as 

AEr,,{t) + (2-4) 

Here, is the relaxation functioiJ^ defined as ds ^m(s) in terms of the 

response function of Umif), = ((V^)['^m(t),^^(0)])^^. The relaxation function is 

independent of temperature, and the value of AEm{0) — AEm{oo) gives the Stokes shift 
magnitude.^ Under the assumption of Gaussian fluctuation and harmonic vibration, the 
Stokes shift is given as double the total reorganization energy and, hence, the relaxation 
function satisfies 


4/m(0) — 2(hAm,env + hAm^vib), (2-5) 

where hXm,env and hAm,vib denote the environmental and vibrational reorganization energies 
of the m-th pigment, respectively. The spectral density, is defined as the imaginary 

component of the susceptibilitj^ and, hence, it is expressed in terms of as 

poo 

.Jmioj) = OJ / dt'^mit) COS Ut. (2.6) 

Jo 

The fluctuations in the electronic transition energies are described by the symmetrized cor¬ 
relation function of iimit), defined as Emit) = ((l/2){-Um(t), ■Um(0)}+),„g. In the quantum 
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mechanical treatment, the Fourier transform of the symmetrized correlation function is ex¬ 
pressed with that of the relaxation function,^ specihcally, = {huj/2) coth(/?ha;/ 2 )d'mM 

In the classical limit of coth(/5ha;/2) ~ 2k-QT/fyjj, therefore, the correlation function of the 
fluctuations can be expressed in terms of ^mif) as 


{Um{t)Ura{^))c\ = k^T ■ 


(2.7) 


Equations (2.6) and (2.7) suggest that the relaxation function, may be the key compo¬ 

nent for characterizing the relevant environmental dynamics and intramolecular vibrational 
motion. 

The relaxation function can be separated into two components: an overdamped part 
originating from the environmental reorganization and an underdamped part induced by 
the intramolecular vibrations, i.e.. 


^m(t) = ^m,env(t) + ^m,vib(t), (2.8) 

with Tm,env( 0 ) = 2hXm,eiw and Tm,vib( 0 ) = 2 hAm„vib- In order to focus on the timescale 
of the environmental dynamics affecting the electronic transition energies, we model the 
environmental component as an exponential decay form, with 




m,env 


(t) 2 fiAenv 6 


■Tenv^ 


(2.9) 


where 7 corresponds to the timescale of the environmental reorganization dynamics. Here, 


it should be noted that the initial behavior of Eq. (2.9) is coarse grained, as was commented 
in Ref. [571 Specihcally, the value of dt'^rn,eiw(t)\t=o = —*hm,env(0) should vanish by dehni- 
tion. However, as long as hnite timescales of the environmental dynamics are discussed, 
this coarse-grained nature of the exponential decay form does not cause fatal defects. For 
simplicity, we consider a single intramolecular vibration on each of the pigments, with fre¬ 
quency, (Uvib, and the Huang-Rhys factor, S. We model the relaxation function originating 
from the vibration, using the Brownian oscillator modeP^ with the vibrational 

relaxation rate, 7 vib, such that 


9 'm,vib(b ^ 2/iA„ib ( cosWvib^ + siliWvibi) e 


a^vib / 

where cDvib = (a^vib ~7vib)^^^ ^-^vib = hcUvib^ have been introduced. Equation (2.6) gives 
the corresponding spectral density, which is reduced to 


7m,vib(ca) = 7ihSu^[6{u - Wvib) - (5(a; -t- cjvib)], 


( 2 . 11 ) 
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in the slow vibrational relaxation limit, i.e., 7vib —t 0. 

An adequate description of the EET dynamics is given by the reduced density operator, 
Pexc(^), that is, the partial trace of the total PPG density operator, pppc(t), over the envi¬ 
ronmental and nuclear DOFs: Pexc(t) = TrenvPppc(t)• The Gaussian nature of the electronic 
excitation energy fluctuations, Um, allows a formally exact equation of motion to be obtained, 
that can describe the EET dynamics under the influence of the environmentally induced fluc¬ 


tuation and dissipation. The nature of these factors is described by Eq. (2.8)-(2.10). Hence, 
the EET dynamics under the influence of the environmental/nuclear dynamics and the cor¬ 
responding 2D electronic spectra can be described in a numerically accurate manner. The 
technical details of this calculation are given in Appendix and Appendix In general, 
vibronic exciton states are introduced as energy eigenstates of the vibronic hamiltonian, 
which are independent of any environmental factors such as temperature and fluctuations. 
However, the correlation function approach in this work does not treat mixed electronic- 
vibrational states as the energy eigenstates, allowing us to describe the influence of environ¬ 
mentally induced fluctuations on the electronic-vibrational quantum mixture. The validity 
of employing the correlation functions to describe the quantum superposition between the 
electronic and vibronic transitions is demonstrated in Appendix 


III. RESULTS AND DISCUSSION 

In this section, we present and discuss the numerical results, in order to quantify the 
impact of the environment upon the quantum mechanically mixed electronic and vibrational 
states on the EET dynamics and 2D electronic spectra of a model dimer. We focus on a dimer 
which produces little beating of electronic origin in the absence of vibronic contributions. 

The studied dimer is inspired by BChl 3 and 4 (pigments 1 and 2, respectively) in the 
FMO protein of Chlorobaculum tepidum^^^^ which serve as the two lowest energy exci¬ 
ton states in the single-excitation manifold. The low-frequency vibrational modes of BChl 
in solutions or protein environments have been investigated in a number of experiments, 
demonstrating that the strongest vibrational mode is found at approximately 180 cm“^.S2ISl 
Thus, to describe the effects of the Franck-Condon active vibrational modes, we consider 
a single vibrational mode with frequency Wvib = 180 cm“^ in this work. For the sake of 
simplicity, the Franck-Condon transition energy of each pigment and their electronic con- 








pling are set to f2i = 12350 cm“\ f22 = — 150 cm“^, and J12 = —50cm“^, respec¬ 

tively. The transition dipole moment directions are set to 9 12 = 90°. In this situation, 
the gap between the two electronic energy eigenstates resonates with the vibrational fre¬ 
quency, [(hli — 0 . 2 )“^ -1-47^2]^^^ — '^vib, and hence, the effects of the vibrational mode are 
expected to be maximized under the given conditions.^ In general, the Huang-Rhys factors 
of BChl are small,!^ and the factor associated with the vibrational mode in this work is 
set to S' = 0.025,^ which is within the range of the experimentally measured values. ^42113] 
The protein-environmental reorganization energy and reorganization time constant are set 
to Aenv = 35cm“^ and 7^^ = 100 fs,[I3lMil ^^d the vibrational relaxation rate is 7^,|, = 2ps. 

In what follows, we let |x^) denote the n-th vibrational level of the a-th electronic state 
in the m-th pigment. At low temperatures satisfying > k-^T, the system can be assumed 
in the \^ig)\x\o) ® l9^2g)|Xgo) state before photoexcitation and, therefore, the two vibronic 
transitions, |95ig)|Xgo) ^ Wie)\xlo) and \^ 2 g)\xlo) ^ \^ 2 e)\xli), can resonate and quantum 
mechanically mix. However, one should not overlook the fact that other resonances may 
be induced as time progresses after the photoexcitation. For example, once the \^ie)\x\o) 
state is created, it is possible for resonance to occur between \^ig)\x]ji) ^ l‘^ie)|Xeo) and 
\V 2 g)\x^go) ^ lv^2e)|Xeo)- ^or the sake of simplicity, we consider the vibrational mode in 
pigment 2 only, as shown in Fig. and we focus on the resonance between the electronic 
transition, \(pig) O \^ie)i of pigment 1 and the vibronic transition, \^ 2 g)\x^go) ^ |v?2e)|xD; 
of pigment 2, which is induced by the coupling ^i2(XeilXgo) = —T12exp(—5/2)^/^. 

A. Population dynamics 

In this subsection, we explore the influence of intramolecular vibrations upon the spatial 
dynamics of electronic excitations in photosynthetic EFT. 

In order to clarify the presence of the vibronic resonance and its influence upon the EET 
dynamics. Fig. presents the time evolution of the energy donor (pigment 1) population 
neglecting environmental factors at temperatures of (a) 77 K and (b) 300 K. In each panel, 
the red line represents the dynamics under the influence of the vibrations, whereas the green 
line gives the dynamics without the effects of the vibrations (i.e., the unitary time evolution) 
as a reference. Both panels indicate that intramolecular vibration dramatically affects the 
EET dynamics in the absence of the environmentally induced fluctuations and dissipation. 
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despite the small Huang-Rhys factor of S' = 0.025. The population dynamics affected by 
the vibrations involve two oscillating components: faster oscillation with small amplitude 
and slower oscillation with large amplitude. By comparing the red and green lines, we can 
observe that the faster oscillating component indicates coherent EET between the electronic 
transition, \(pig) -f-)- \^ie)i and the 0 — 0 vibronic transition, \^2g)\')dgo) \^2e)\xlio)i which is 
induced by the coupling constant, ^i2(XeolXgo) — >^12 exp(—S'/ 2 ) = 0.987Ji2 for S = 0.025. 
On the other hand, the slower oscillating component arises from the interaction between the 
electronic transition, \^ig) -H- |(pie), and the 0 — 1 vibronic transition, |(p 2 c/)|Xgo) ^ \V2e)\')dei)) 
which is induced by the coupling constant, <7i2(XeilXgo) ~ — J 12 exp(—S'/2)\/5 = —0.156J 12 . 
Moreover, the large amplitude of this slower component indicates a large mixing angle, which 
is the resonance between the electronic and vibronic transitions. 


Figure 1^ shows the time evolutions of the energy donor (pigment 1) population affected 
by the intramolecular vibrations under the influence of the environmentally induced fluc¬ 
tuations and dissipation at temperatures of (a) 77 K and (b) 300 K. Similar to Fig. the 
red lines represent the dynamics in the presence of the vibrations, whereas the green lines 
show the dynamics without the vibrations as a reference. The EET dynamics in Fig. 
show no wavelike motion in contrast with those in Fig. However, it should be noted that 
the dynamics involve two distinct mechanisms. Comparing Fig. [3] with Fig. 1^ one recog¬ 
nizes that the short-term behaviors {t < 7 ^^ = 100 fs) arise from electronically coherent 
motion between the electronic transition of pigment 1 and the 0-0 vibronic transition 
of pigment 2 , and the subsequent hopping dynamics follow. As demonstrated in Ref. [57], 
the electronic excitation is delocalized just after the photoexcitation even in the incoherent 
hopping regime. As time progresses, the dissipation of reorganization energy proceeds and 
the excitation becomes localized. The vibrational modes contribute to acceleration of the 
EET, albeit only slightly, because the vibronic resonance adds a new EET channel in com¬ 
parison with the case without the vibrational mode. However, the overall behaviors of the 
EET dynamics are dominated by the environment, and the contribution of the vibrational 


modes decreases with increasing temperature. As can be seen in Eqs. (2.7) and (2.9), the 
root mean squared amplitude of the environmentally induced fluctuations in the electronic 
energy is given as ( 2 hAenv/^BT)^/^, which is evaluated as 61.2 and 120.8 cm“^ at temperatures 
of 77 and 300 K, respectively. On the other hand, the inter-site coupling constant inducing 
the vibronic resonance is evaluated as 7i2(XeilXgo) = 7.80cm“^ for S = 0.025. Hence, the 
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inter-site coupling causing the vibronic resonance is one order of magnitude smaller than 
the fluctuation amplitudes and, therefore, the vibronic resonance does not play a signihcant 
role in the EET dynamics under the influence of the environmentally induced fluctuations. 


B. Beating in two-dimensional electronic spectra 

In this subsection, we discuss 2D electronic spectra, in order to explore the relation¬ 
ships between the spatial dynamics of the electronic excitations presented in the preceding 
subsection and the information content of the 2D electronic spectra. 

Figures]^ and|^ show calculated nonrephasing 2D spectra at 77K without and with the 
intramolecular vibration, respectively. These calculations were performed using the same 
equation of motion and parameters as in Fig. The lower-left and upper-right diagonal 
peaks are labeled DPI and DP2, respectively. The amplitudes of the diagonal cuts are 
also presented as a function of the waiting time, ^ 2 , in Fig. and|^. Although the EET 
dynamics in Fig.|^and the 2D spectra at ^2 = 0 do not exhibit a signihcant change originating 
from the vibrational mode, Figs.|^ and|^ exhibit a qualitative difference. Figure]^ shows 
monotonous decays of the diagonal peaks, whereas Fig. exhibits peak beating persisting 
for up to a minimum of 2ps. Apparently, this beating is caused by the presence of the 
vibrational mode. 

To clarify the origins of the beatings in Fig. ii. we decompose the nonrephasing 2D 
spectra to the stimulated emission (SE), ground state bleaching (GSB), and excited state 
absorption (ESA) contributions. Figure shows the time-evolutions of the diagonal cuts of 
the individual contributions. The time-evolutions of the amplitudes of the diagonal peaks, 
DPI (wi = uis = 12,166cm“^) and DP2 (cui = = 12,350cm“^) in the contributions 

are also shown in Fig. The SE and ESA contributions exhibit beating corresponding 
to quantum superposition between the vibronic exciton states, |eo) and |ei*"), in Fig. 
whereas the beating in the GSB contribution arises from quantum superposition between 
the vibrational states in pigment 2, l^fg) and As shown in Figs. I and 1^ the SE, 

GSB, and ESA pathways exhibit similar contribution levels to the diagonal peak beatings 
in the nonrephasing 2D spectra, although the SE and ESA pathways exhibit larger beating 
amplitudes in the short-term region (~ 300 fs), which originate from the electronic coherence 
in the single-excitation manifold. The present observation corroborates the recent theoretical 
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argument by Plenio, Almeida, and Huelga,!^ which suggested that the SE, GSB, and ESA 
contributions in the rephasing 2D spectra could be of the same order in realistic parameter 
ranges with the use of a simplihed Markovian quantum master equation. 

In order to clarify the extent to which the quantum mixture between the electronic tran¬ 
sition of pigment 1 and the vibronic transitions of pigment 2 contributes to the observed 
beating in the nonrephasing 2D spectra, we explore the effects of the environmentally in¬ 
duced fluctuations on beating amplitude. For this purpose, we address beating in the GSB 
pathways, that is, the vibrational coherence between Iqq) and \gi) in the electronic ground 
state. This is because the GSB pathway is unaffected by the energy transfer processes 
in the single-excitation manifold and, thus, it allows one to investigate interplays between 
the quantum mixture and the environmentally induced fluctuations. The beating ampli¬ 
tudes in the GSB contribution observed in the nonrephasing 2D spectra are determined 
as being the difference between the hrst and second extremal values after G = 500 fs, as 
indicated by the dashed-line box in Fig. |^, where the environmental reorganization dynam¬ 
ics have already been completed and the beating characteristics can be extracted cleanly. 
Figure shows the beating amplitude along the diagonal line of the GSB contribution 
for various values of cUvib- From the plots, we hnd two characteristic beating behaviors in 
the GSB contribution of the nonrephasing 2D spectra. Firstly, a local maximum of the 
beating amplitude is observed at the DPI position (cui = 0 J 3 = 12,166cm“^), indepen¬ 
dently of the vibrational frequency. This beating is caused by the ground state vibrational 
coherence created via the lower exciton, |eo), that is, through the Liouville pathway of 
l5'o)(5'o| -t |eo)(fi'o| -t |5'i)(5'o| -t |eo)(5'o| -t |fi'o)(fi'o|- Secondly, another local maximum is 
found at various locations in the vicinity of DP2 (cui = U 3 = 12,350cm“^), depending on 
the vibrational frequency. This beating originates from the ground state vibrational coher¬ 
ence created via the upper vibronic excitons, |ei*"), that is, through the Liouville pathways 
of |5'o)(5'o| -t |ei^)(5'o| -t |5'i)(5'o| -t |e^)(5'o| -t |5'o)(5'o|- The transition energy of the vi¬ 
bronic excitons, |ef), differs according to the vibrational frequency and, hence, the location 
of the beating amplitude maxima vary. Moreover, the mixing angle between the electronic 
transition of pigment 1 and the vibronic transition of pigment 2 depends on the vibration 
frequency, which causes variations in the beating amplitudes. This fact leads to the concept 
that the degree of quantum mixing between the electronic transition of pigment 1 and the 
vibronic transition of pigment 2 may be visualized as the beating amplitude of the ground 
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state vibrational coherence. 

Figure plots the beating amplitude of DPI and the local maxima values of the beating 
amplitude in the vicinity of DP2 as a function of Wvib- The beating at the DPI location 
originates primarily from the ground state vibrational coherence created via |eo), which is 
minimally dependent on the quantum mixture between the electronic and vibronic tran¬ 
sitions. However, the beating amplitude of DPI decreases with the increasing vibrational 
frequency. This behavior is explained as being due to the interference between the beat¬ 
ing of DPI and that in the vicinity of DP2, as shown in Fig. On the other hand, the 
beating amplitude near DP2 exhibits a maximum in the vicinity of cjvib = 180 cm“^. These 
hndings reflect the degree of quantum mixing and, thus, the degree of intensity borrowing. 
Tiwari, Peters, and Jona^ demonstrated that the beating amplitude of a cross-peak in the 
rephrasing GSB pathway may be enhanced over a wide vibrational frequency range around 
the vibronic resonance. Their hnding is consistent with those of the present study that 
adequately describes the impact of the environmentally induced fluctuation and dissipation 
upon the quantum mixing. Figure [^implies that the quantum mixing between the vibronic 
transitions in BChl in the FMO protein, which induces vibronic coherence in the electronic 
ground and excited states, can be rather robust. This holds even under the influence of the 
environmentally induced fluctuations, at least at cryogenic temperatures. 

However, the root mean squared amplitude of the environmentally induced fluctuations 
in the electronic energy is estimated as (2hAenvfcBT)^/^ using Eqs. (2.7) and (2.9), and the 
quantum mixing between the strong electronic and weak vibronic transitions may be more 
fragile against the protein-induced fluctuations at higher temperatures. Figure shows the 
beating amplitude of DPI and the local maximum value of the beating amplitude in the 
vicinity of DP2 as a function of temperature for Wvib = 180cm“^. The beating of DPI 
originates from the ground state vibrational coherence created via |eo), independently of 
the quantum mixture and, thus, it exhibits minimal temperature dependence as well as the 
vibrational frequency dependence apparent in Fig. However, the beating amplitude in the 
DP2 region is strongly dependent on temperature. This indicates that the environmentally 
induced fluctuations of larger amplitude at higher temperature eradicate the quantum mixing 
between the strong electronic transition of pigment 1 and the weak vibronic transition of 
pigment 2. At the high temperature limit, the quantum mixture and, thus, the intensity 
borrowing mechanism, diminish. Hence, the beating amplitudes of DPI and DP2 converge 
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towards each other. Figure shows the time-evolutions of the amplitudes of the diagonal 
cuts of nonrephasing 2D electronic spectra at 300 K for (a) a coupled dimer without the 
intramolecular vibration in pigment 2 and (b) a coupled dimer with the intramolecular 
vibration. The intensity borrowing mechanism is eradicated at 300 K, as indicated in Fig. 
and the vibrational contribution in Fig. pTj p almost vanishes as a result. 

The fluorescence line-narrowing measurement of the FMO protein reveals that a vibra¬ 
tional mode at 117cm“^ and several modes between 167 — 202 cm“^ have relatively strong 
Huang-Rhys factor magnitudes.S^l As shown in Fig. the beating amplitude in the GSB 
contribution is enhanced over a wide range of vibrational frequencies near cUvib = 180 cm“^ 
at cryogenic temperature, i.e., 77 K. This behavior supports the recent theoretical result 
reported by Tempelaar, Jansen, and Knoester,^ who showed that the vibrational modes of 
117 cm“^ and 185 cm“^ induce long-lasting beating of the exciton 1 — 3 cross peak in the 
rephasing 2D spectra of the FMO protein. This was also investigated by Panitchayangkoon 
et a/.® However, the amplitudes of the beats induced by the electronic-vibrational quantum 
mixture decrease with increasing temperature, as shown in Figs. Moreover, the vibra¬ 
tional contribution almost vanishes, as shown in Fig. 10 Panitchayangkoon et alW showed 
that quantum beats in the FMO complex persist for a minimum of 1.5 ps and 300 fs at 77 
and 277 K, respectively, whereas Ishizaki and Flemin^^® theoretically predicted that elec¬ 
tronic coherence in the FMO complex could persist for 700 and 300 fs at 77 and 300 K, 
respectively. The experimental data showed good agreement with the theoretical prediction 
regarding the electronic coherence at physiological temperature; however, the experimental 
data exhibited a signihcantly longer quantum beat lifetime in comparison with the theoret¬ 
ical result. Hence, the temperature dependence presented in Figs. ii and indicates 
that the long-term behavior of the experimentally observed long-lasting quantum beats is 
primarily vibrational in nature. 


IV. CONCLUDING REMARKS 

In this work, we have explored the impact of environmentally induced fluctuations upon 
the quantum mechanically mixed electronic and vibrational states of pigments in terms of 
EFT dynamics and the corresponding 2D electronic spectra. For this purpose, we have 
addressed a model dimer extracted from the FMO protein and have performed numerical 
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calculations of EET dynamics and the corresponding 2D electronic spectra in an accurate 
manner. 

At cryogenic temperature, the resonance between the strong 0 — 0 and weak 0 — 1 vibronic 
pigment transitions are rather robust, even under the influence of protein-induced fluctua¬ 
tions, despite the use of a small Huang-Rhys factor, i.e., S = 0.025. This resonance results 
in quantum beats of vibrational origin in both the electronic ground and excited states in 
the 2D electronic spectra via the intensity borrowing mechanism. This is not specihc to 
the hne resonance between the vibronic transitions, and the beating originating from the 
vibrational coherence can be enhanced over a rather wide range of vibrational frequencies 
around the resonance condition, at cryogenic temperature at least. However, the amplitudes 
of the environmentally induced fluctuations in the electronic excitation energy increase with 
increasing temperature. The present calculations demonstrate that the quantum mixing 
between the strong 0-0 and 0 — 1 vibronic transitions is more fragile against fluctuations 
at higher temperatures, and the quantum beats originating from the vibrational coherence 
in 2D electronic spectra almost vanish at physiological temperatures. 

Furthermore, we have investigated EET dynamics using the same dynamic equation and 
parameters as in the calculations for the 2D electronic spectra. In the absence of the en¬ 
vironmental effects, the resonance between the vibronic transitions dramatically alters the 
EET dynamics, despite the small S (= 0.025). However, under the influence of the envi¬ 
ronmentally induced fluctuations and dissipation, the vibronic resonance contributes only 
slightly to acceleration of the electronic energy transfer, both at cryogenic and physiologi¬ 
cal temperatures, although the 2D electronic spectra clearly exhibit quantum beating that 
is enhanced by the vibronic resonance at cryogenic temperature. One of the remarkable 
benehts of nonlinear multidimensional spectroscopic techniques is their sensitivity to details 
of the systems and dynamics under investigation. The calculated 2D electronic spectra in 
the present study indicate vanishingly small effects of the vibronic quantum mixing upon 
the EET dynamics and present them as long-lived quantum beating. In other words, this 
suggests that the concept that oscillatory behavior in 2D data provides genuinely relevant 
information on the systems and dynamics under investigation should be carefully considered. 
Recently, Fuller et a/.® and Romero et revealed the presence of long-lived quantum beats 
in the photosystem H reaction center through the use of 2D electronic spectroscopy. Both 
works demonstrate that the observed long-lived beats are induced by the electronic and 
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vibrational resonance between the electronic exciton and the primary charge transfer states, 
and they suggest that the electronic-vibrational resonance may represent an important de¬ 
sign principle for enabling the high quantum efficiency of charge separation in the reaction 
center. These findings are indeed intriguing. However, the reorganization energies and, thus, 
the protein-induced fluctuations associated with charge transfer states are generally large 
in comparison to those for the electronic exciton states and, hence, further investigation of 
the influence of the electronic-vibrational resonance upon the charge transfer dynamics is 
required. This will be the subject of future study. 
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Appendix A: Equation of motion 


In order to derive the equation of motion to describe EET dynamics, we organize the 
product states in order of elementary excitation number. The overall ground state with 
zero excitation reads \g) = \(pig)\g:> 2 g)■ The presence of a single excitation at pigment 1 is 
expressed as |1) = \Lpie)\ip 2 g), whereas the other is |2) = \Lpig)\Lp 2 e)■ Since the intensity of 
sunlight is weak, the single-excitation manifold is of primary importance under physiological 
conditions. However, nonlinear spectroscopic techniques such as 2D electronic spectroscopy 
can populate some higher excitation manifolds, e.g., the double-excitation manifold com¬ 


prising 112) = \ipif,)\if 2 e)■ Hence, the PPG Hamiltonian in Eq. (2.1) can be recast as 


H 


PPG 


Hexc T 77e 


Hexc—envi 


01 ) 
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with 


Hexc Eg 15 ^) (^ 71 

2 2 

h-t/ mn\fE) (jl 

m=l m=l n^m 


+ [Eg + hrii + 112) (12|, 

@2) 

2 

Henv ^ ^ HyYigi^Xm) 

m=l 

@3) 

2 

Hexc-env = Um\m) {m\ + (ui M 2 ) 112) (12|. 

@4) 


m=l 


In the above, the excitation-vacuum energy, Eg = J2m{^mg{xm))mg, has been introduced; 
however, we set Eg = 0 without loss of generality. In order to derive the reduced density 
operator to describe the EET dynamics, we suppose that the total PPG system at the 
initial time, t = 0, is in the factorized product state of the form, Ppp^ = Pexc Hm This 
factorized initial condition is generally considered to be unphysical in the literature on open 
quantum systems,!^ since it neglects the inherent correlation between a system of interest 
and its environment. In electronic excitation processes, however, this initial condition is of no 
consequence, because it corresponds to the electronic ground state or an electronic excited 
state generated in accordance with the vertical Franck-Condon transition. The reduced 
density operator for the electronic excitation, Pexcit), is given by Pexcit) = Trenv[G(f)ppp(^], 
where G{t) is the retarded propagator of the total PPG system in the Liouville space. The 
Gaussian property of Um enables one to derive a formally exact equation of motion a^ 

^Pexc(i) = / dsK^^\t, s)Pexcit) + T+ f dsK[l\t, s)PexciG, 

m=l “^0 



where the interaction representation has been employed with respect to Hexc + He 
integration kernel for the single-excitation manifold, Km\t,s), is given by 




Dmit - s)Kn(s)'^ - i^^mit “ 


The 


06 ) 


where Vm = |m)(m| has been introduced, and Emit) and are the symmetrized cor¬ 

relation function and the response function of the collective energy gap coordinate u{t), re¬ 
spectively. On the other hand, the integration kernel used to describe the double-excitation 
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manifold, Km (t, s), is given by 




fll2(i - s)Vi2(s)>‘ - 2^4 i 2(( - s)n2(s)‘ 


07 ) 


with Duit) = Di(t) + D 2 {t) and $i 2 (t) = *hi(t) + <h 2 (^)- The superoperator notation has 
been introduced, 0^02 = O 1 O 2 — O 2 O 1 and 0^02 = O 1 O 2 + O 2 O 1 , for any operators Oi and 


O 2 . In Eq. (A 5), the chronological time ordering operator, T+, resequences and mixes the 
superoperators, Vm(i)^ and Vm{t)°, comprised in Km\t,s) and Pexc{t) and, hence, Eq. (A 5) 


is difficult to solve. In general, Eq. (A 5) can be recast exactly or approximately into the 
form, 

^Pexc(t) = y"T+ @8) 

at ^ Jq 

where and Ca are arbitrary operators and is an arbitrary complex number, and 


do not necessarily have physical meaning. In this situation, Eq. (A 5) is recast into the 
multidimensional recurrence form a^SZESEZl 
d 


—a{ni,n 2 , • • •) = -{t^exc + riibi + 71262 + • • • )d{ni,n 2 , ...) 

+ Aid{ni + 1, 772, . . . ) + ^2d(77i, 772 + 1, ...) + .. . 

+ niCid{ni - 1,772, . . . ) + 772C'2d(77i,772 - 1, , 


09 ) 


where £exc is the Liouvillian corresponding to the Hamiltonian, H^xc in Eq- (A 1). In the 


above recurrence formula, only the d(0, 0,...) element is identical to the reduced density 
operator, pexc{t), while the other terms are auxiliary operators which do not necessarily 


have physical meaning. The recurrence in Eq. (A 9) continues to infinity. However, the 


numerical calculations can converge at a finite stage. To terminate Eq. (A 9) safely, we 


replace Eq. (A 9) with 


— (j(77i, 772, . . . ) = -7£excd(771,772, . . . ), 


010 ) 


for the integers ( 771 , 772 ,...) satisfying S> Uc/ minQ,(6a), where Uc is a characteristic 

frequency for 

Appendix B: Nonrephasing third-order response function 

The nonrephasing 2D spectrum is calculated as 


*S'nr(i^3 )^2, 1 ^ 1 ) — Re / ^2,H), 

Jo Jo 


01) 
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where the nonrephasing response function, -RNR(t 3 , ^ 2 , ti), is expressed 

-Rnr(^ 3 ,^ 2 ,^i) = Tr ^il^G(t 3 )-fi^G{t 2 )-fi^G{ti)-fi^p^ppfj'^ . @2) 

The operators, /t_j, and /t^, are defined as 

h-s> = /^i|l)(5'l + /^2|2)((y'| + /ii|12)(2| + /i2|12)(l|, @3) 

/i^ =/ii|^)(l| +/i2|(7)(2| +/ii|2)(12| +/i2|l)(12|, @4) 


where is the transition dipole of the m-th pigment. We assume the Condon approxima¬ 
tion and, thus, is independent of the environmental and nuclear DOFs. For simplicity, 
we set fim = 1 and /xi ■ = cos 6 * 12 , with 6*12 being the angle between the two transition 


dipoles. Note that the directions of the subscript arrows in Eq. (B 2) correspond to those 
of the arrows in the double-sided Feynman diagrams depicting the nonrephasing response 


functions. Regarding the means of calculating the response function in Eq. (B2) using 
the equation of motion approach given in Appendix we refer to Sec. 3 in Ref. |6HI and 
Sec. 5.5 in Ref. Eni We note that numerical results of the third-order response function in 


Eq. (B 2) do not depend on employed representations (e.g. the site representation or the 


energy eigenstate representation) because Eq. (B 2) is expressed as the trace of operators. 


Appendix C: Validity of the correlation function approach to describe vibronic 
excitons 


Here, we comment on the validity of employing the correlation functions and the cor¬ 
responding spectral density to describe the mixed electronic-vibrational state, instead of 
employing the vibronic exciton model.l^ F7 | 28 | 46 | consider the expression of Forster’s EET 
rat d^°*^^ * 


'^A-s-D ~ "^AD 




( 01 ) 


where /a[ a)] and EdIu] are the absorption line-shape of the energy acceptor molecule (A) 
and the emission line-shape of the donor (D), respectively. For simplicity, we employ the 


Jm, given in Eq. ( 2 . 11 ) and assume the zero-temperature limit. In this situation, 
Eq. (|gg is recast into 


du 




v=0 v'=0 


102 ) 
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where 


i(x?„i\'r„>p = (03) 

= 7t6{u - {Qa - •S'cuvib + vtJvib)), 04) 

E^ \uj] = 7 t6{u - (r^D - -S'wvib - ^^'c<;vib))• 05) 


In the above, is the Franck-Condon factor associated with the vibronic transi¬ 

tion from the n-th vibrational level of the a-th electronic state, \x^v)i to the n'-th vibrational 
level of the a'-th electronic state, \x^iv')i of the m-th pigment. Equation (jC 2) indicates that, 
in the case of hlA + '^i^vib = the two vibronic transitions, |(^De)|x°o) \V'Dg)\x^v') 

and \^Ag)\x%) WAe)\xtv)^ resonate, and excitation transfer then occurs between them 
with the coupling constant 


JATi{Xev\X%){X^v'\X^o)- 06) 

This is consistent with the approach using the vibronic exciton Hamiltonian mo del, 
and, thus, the present correlation function approach certainly describes the interaction be¬ 
tween the vibronic transitions. Furthermore, it should be emphasized that the current 
approach is capable of describing vibrational relaxations, vibrational states in both the elec¬ 
tronic ground and excited states, and the influence of the environmental dynamics in a 
natural and consistent manner. 
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FIG. 1. (a) Illustration of quantum mechanical mixing among vibronic transitions of pigments 1 
and 2 in a weakly coupled hetero-dimer. Arrows are used to depict strong vibronic transitions 
only. For simplicity, only a single vibrational mode is considered in pigment 2. The exciton state, 

I Co), originates primarily from the vibronic state of pigment 2, |(^2e)|Xeo)- Resonance or strong 
quantum mixing between the electronic state of pigment 1, \^ie), and the vibronic state of pigment 
2, \if 2 e)\'X^i) 1 creates the two vibronic excitons, |ef), in the absence of the environmentally induced 
fluctuation and dissipation. In the vibronic exciton picture, the overall ground state, l^o); is given 
by Iv’is)® l¥^2g)|Xgo)) while |( 7 i) is identical to |<^2g)|Xgi)- (b) The relaxation functions in Eqs. (2.8)- 
( |2.10 ) and the corresponding spectral densities, Eq. (2.6), employed for modeling the environmental 
dynamics and intramolecular vibrational motion affecting the electronic transitions of pigments 1 
and 2. 
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Without the environment 


(a) temperature 77K (b) temperature 300K 



time (fs) time (fs) 


FIG. 2. Time evolution of the energy donor population affected by intramolecular vibration of 
pigment 2 (acceptor) in the absence of protein environment effects at temperatures of 77 and 300 K. 
In each panel, the red line represents the dynamics under the influence of the vibration, while the 
green line shows the dynamics without the vibration effects, i.e., the unitary time evolution. The 
vibrational frequency, the vibrational relaxation time constant, and the Huang-Rhys factor are set 
to cuvib = 180cm“^, = 2ps, and S = 0.025, respectively. 


With the environment 


(a) temperature 77K (b) temperature 300K 



time (fs) time (fs) 


FIG. 3. Time evolution of the energy donor population affected by intramolecular vibration of 
pigment 2 (acceptor) in the presence of environmental effects at temperatures of (a) 77 and (b) 
300 K. In each panel, the red line represents the dynamics under the influence of the vibration, 
whereas the green line shows the dynamics without the vibration effects. The environmental 
reorganization energy and the reorganization time constant are set to Aenv = 35cm“^ and 7 ^^ = 
100 fs, respectively. The other parameters are the same as in Fig. 
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FIG. 4. Nonrephasing 2D spectra at 77 K of (a) a coupled-dimer without intramolecular vibration 
and (b) a coupled-dimer with intramolecular vibration. The waiting time is t 2 = Ofs. The nor¬ 
malization of contour plots (a) and (b) is such that the maximum value of each spectrum is unity, 
and equally spaced contour levels (0, ±0.02, ±0.04,...,) are drawn. The lower-left and upper-right 
diagonal peaks are labeled DPI and DP2, respectively. Panels (c) and (d) show the amplitudes of 
the diagonal cuts of 2D spectra (a) and (b) as functions of the waiting time, t 2 -> respectively. The 
calculations were performed with the parameters employed in Fig. 
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FIG. 5. Time evolutions of the diagonal cuts of (a) the stimulated emission, (b) the ground state 
bleaching, and (c) the excited state absorption contributions to the nonrephasing 2D spectra given 
in Figs. & and P- The normalization of the contour plots is such that the maximum value of 
Fig.p is unity, and equally spaced contour levels (0, ±0.02, ±0.04,...,) are drawn. 


27 























waiting time (fs) 


FIG. 6 . Time evolutions of the diagonal peak amplitudes, DPI (cui = L03 = 12,166cm“^) and 
DP2 (wi = W 3 = 12,350cm“^), of (a) the nonrephasing 2D spectrum, (b) the stimulated emission 
contribution, (c) the ground state bleaching contribution, and (d) the excited state absorption 
contribution. The normalization of the plots is such that the maximum value of Fig. & is unity. 
In order to clarify the relationship between the beating amplitude magnitudes in the individual 
contributions, the peak amplitudes are shifted. 
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FIG. 7. Beating amplitude distribution along the diagonal line of the GSB contribution in the 
nonrephasing 2D spectra for various values of the vibrational frequency, Wvib- The temperature is 
77 K. The normalization of the individual plots is such that the maximum value of Fig. Ep is unity. 



FIG. 8 . Beating amplitude of DPI (wi = 013 = 12,166cm“^) and the local maximum value of the 
beating amplitude in the vicinity of DP2 (wi = 003 = 12,350 cm“^) as a function of the vibrational 
frequency, cuvib- The temperature is 77 K. The normalization of the individual plots is such that 
the maximum value of Fig. is unity. 


29 









0.14 



FIG. 9. Beating amplitude of DPI (cui = io^ = 12,166cm“^) and the local maximum value of the 
beating amplitude in the vicinity of DP2 (wi = wa = 12,350cm“^) as a function of temperature. 
The vibrational frequency is fixed to cuvib = 180 cm“^. The normalization of the individual plots 
is such that the maximum value of Fig. is unity. 
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FIG. 10. Time evolutions of the diagonal cut of the nonrephasing 2D spectra at 300 K of (a) 
a coupled-dimer without intramolecular vibration and (b) a coupled-dimer with intramolecular 
vibration. The normalization of contour plots (a) and (b) is such that the maximum values of 
Fig.|4}r and are unity. Equally spaced contour levels (0, ±0.02, ±0.04,...,) are drawn. 
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